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We report on studies of quantum turbulence with second-sound in superfluid ^He in which the 
turbulence is generated by the flow of the superfluid component through a wide square channel, the 
ends of which are plugged with sintered silver super leaks, the flow being generated by compression 
of a bellows. The super leaks ensure that there is no net flow of the normal fluid. In an earlier 
paper [Phys. Rev. B, 86, 134515 (2012)] we have shown that steady flow of this kind generates 
a density of vortex lines that is essentially identical with that generated by thermal counterflow, 
when the average relative velocity between the two fluids is the same. In this paper we report 
on studies of the temporal decay of the vortex-line density, observed when the bellows is stopped, 
and we compare the results with those obtained from the temporal decay of thermal counterflow 
re-measured in the same channel and under the same conditions. In both cases there is an initial 
fast decay which, for low enough initial line density approaches for a short time the form t~^ 
characteristic of the decay of a random vortex tangle. This is followed at late times by a slower 
decay, characteristic of the decay of large “quasi-classical” eddies. However, in the range of 
investigated parameters, we observe always in the case of thermal counterflow, and only in a few 
cases of high steady-state velocity in superflow, an intermediate regime in which the decay either 
does not proceed monotonically with time or passes through a point of inflexion. This difference, 
established firmly by our experiments, might represent one essential ingredient for the full theoretical 
understanding of counterflow turbulence. 

PACS numbers: 67.25.dg, 67.25.dk, 67.25.dm 


I. INTRODUCTION 

In this paper we report on an experimental investiga¬ 
tion of the decay of turbulence in a quantum fluid, su¬ 
perfluid ^He, displaying the two-fluid behavior. The nec¬ 
essary introduction to superfluidity and quantum turbu¬ 
lence is given in the next section, together with a review 
on counterflow turbulence, a unique form of turbulence 
existing in superfluid ^He, of which the superflow treated 
in this article represents a special case. The review sets 
this work into a detailed context. Here we summarize the 
general motivation for this project, which is two-fold. 

(i) We wish to continue the study of temporal decay of 
turbulence in quantum fluids, in the spirit of exploring 
similarities and differences with the decay of turbulence 
in classical viscous fluids, as advocated again as an im¬ 
portant goal for the community in the latest review on 
the state of the field [l|. The decay of turbulence is in¬ 
deed a cornerstone in classical turbulence studies because 
the rate of decay (of kinetic energy or vorticity) is related 
to the distribution of energy over the scales of the sys¬ 
tem, which constitutes a key description of a turbulent 
flow [H . In quantum turbulence our experiments can ac¬ 
curately measure the decay of the total length of quan¬ 
tized vortices per unit volume, a well defined quantity 
which, if the detailed spatial distribution of the vortex 
lines is known, can be related to classical vorticity. 

(ii) Our second and more specific motivation is to con¬ 
tinue our investigation of the mechanically-driven turbu¬ 
lent flow of the superfluid component of ^He, of which we 
have reported the steady-flow characteristics in this jour¬ 


nal Q . The study of pure superflow ought to be in princi¬ 
ple simpler than thermal counterflow because the normal 
component is on average at rest. In our previous work 
we have demonstrated that Galilean invariance holds be¬ 
tween steady-state counterflow and superflow turbulence, 
i.e. that, to first order, the turbulence produced when the 
superfluid and normal components of ^He move in oppo¬ 
site direction is the same as when the normal fluid is at 
rest and the superfluid moves past it with the same rela¬ 
tive velocity, despite the fact that the presence of a finite 
channel ought to introduce differences such as a profile 
to the normal component. In this paper we extend this 
work by asking how these two turbulent flows decay in 
time when driving is suddenly switched off. We describe 
emerging similarities as well as differences, which ought 
to deepen our general understanding of the underlying 
physics of quantum turbulence, especially in relation to 
the dynamical state of the normal component. 


II. REVIEW OF COUNTERFLOW 
TURBULENCE 

Quantum turbulence is the turbulence occur¬ 

ring in a superfluid such as the superfluid phases of liquid 
^He and liquid ^He Q- At a finite temperature super¬ 
fluids exhibit two-fluid behavior, a normal viscous fluid 
(composed of thermal excitations) coexisting with an in- 
viscid superfluid component. Flow of the superfluid com¬ 
ponent is strongly influenced by quantum effects, reflect¬ 
ing the origin of superfluidity in Bose condensation. In 
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the case of ^He superfluid flow must be irrotational, ro¬ 
tational motion being possible in a simply-connected vol¬ 
ume only with the formation of topological defects in the 
form of vortex lines, each of which carries a circulation 
oii<i = h/m^lx 10“^ cm^/s, where h is Planck’s con¬ 
stant and m is the mass of a helium atom. Turbulence in 
the superfluid component must therefore in general take 
the form of a complex tangle of vortex lines. A purely 
random tangle involves turbulent energy almost entirely 
on only quantum length scales - scales comparable with, 
or less than, the mean vortex spacing, i = where 

L is vortex line density (length of line per unit volume), 
although local polarization of the vortices can lead to the 
existence of turbulent energy on any larger scale iii- 
Both purely quantum and classical features of turbulence 
can therefore be detected simultaneously in the same 
quantum flow, depending on the length scale at which 
this quantum flow is probed i,®. 

If vortex lines move relative to the normal fluid they 
experience a drag force, called mutual friction Hi- On 
quasi-classical length scales; i.e., scales large compared 
with the superfluid component usually behaves like a 
classical fluid at high Reynolds number. The small kine¬ 
matic viscosity of the normal fluid (of order k ,/6 0 ) 
means that the same is true for the normal fluid on these 
large quasi-classical scales. Thus the two fluids can move 
together with the same velocity fields, mutual friction 
serving only to stabilize this coupled motion. The cou¬ 
pled fluids then behave as a single quasi-classical fluid 
at high Reynolds number. This situation, referred to 
as CO flow ^ obtains quite frequently; for example when 
the fluid is stirred by large scales objects similarly as 
for classical fluids, such as propellers [l 2 |, grids 0 , flow 
through channels [iHi3 or due to various oscillating ob¬ 
jects |lq . 

On quantum length scales (even in coflow) the super- 
fluid motion is strongly influenced by quantum effects. 
The coupled motion therefore cannot be maintained, the 
resulting dissipation due to the mutual friction, combined 
with viscous dissipation in the normal fluid, results in a 
dissipation per unit mass of helium of the form 

e = v' ^ (1) 

where v' is an effective kinematic viscosity of order 
n Bill, [Bill. This type of coupled motion, leading to 
quasi-classical behavior on large length scales, will usu¬ 
ally obtain as long as there is no forced relative motion 
between the two fluids on these scales. We shall refer to 
such forced relative motion as counterflow. It is most eas¬ 
ily imposed with a temperature gradient, the superfluid 
component moving up the gradient and the normal fluid 
moving down, with no net mass flow. This is one special 
case of counterflow, known as thermal counterflow. 

Quite generally, by combining mechanical and thermal 
drive, a rich variety of counterflows, i.e., two-fluid flows 
with different flow ratio of the two components, can be 
generated, representing a very complex superfluid hydro- 
dynamic system. In this paper, we are concerned with 


the special case of counterflow that can be generated by 
forcing helium through a tube, the ends of which are 
closed with superleaks - only the superfluid component 
can flow through the super leaks, so that the average ve¬ 
locity of the normal fluid vanishes. The forced flow is 
conveniently driven by compression of a bellows, and we 
refer to this type of counterflow as bellows-driven super¬ 
flow 0 , or simply superflow. 

Early experiments on thermal counterflow led to 
the idea that in such a flow a self-sustaining random tan¬ 
gle of vortex lines is generated simply by the imposed 
relative motion, the turbulence being essentially homo¬ 
geneous, and a phenomenological equation was derived 
that describes the growth and decay of line density. In 
the steady state the line density is given in terms of the 
relative velocity v = {vs — Vn) by 

L = (2) 

where 7 is a temperature-dependent parameter. Confir¬ 
mation of these ideas, together with an understanding of 
the physical processes involved, came from the pioneer¬ 
ing computer simulations of Schwarz 1^ . which were im¬ 
proved and refined by Adachi et al. 1 ^ . Baggaley and 
Barenghi [25[ and Kondaurova et al. [26| . These simu¬ 
lations were based on the vortex filament model. They 
assumed that the flow of the normal fluid remained lam¬ 
inar, and that the mean flow velocities were spatially ho¬ 
mogeneous and in an unbounded volume. The resulting 
vortex tangle was disordered, so there was no large scale 
turbulent motion in the superfluid component. 

In Ref. 0 we reported the results of a study of the at¬ 
tenuation of second-sound in bellows-driven steady-state 
superfiow. According to the theories that we have de¬ 
scribed so far, the line densities observed in such an ex¬ 
periment should agree with those measured at the same 
value of V as those observed in thermal counterflow. This 
was indeed confirmed to a good approximation in chan¬ 
nels of different cross-sections, with the agreement be¬ 
ing particularly good when counterflow and superflow are 
measured in the same (large) channel 0 . 

There are, however, subtle differences between coun¬ 
terflow and superflow, even in the steady-state. Exper¬ 
iments have shown that there is a critical velocity be¬ 
low which the line density is unmeasurably small, which 
scales with the channel size, D. It was found to be 
roughly temperature independent scaling as D~^l^ for 
superflow, while for counterflow it displays d, D~^ tem¬ 
perature dependent scaling Q. Additionally, another 
critical velocity has been reported by Tough and his co¬ 
workers [23 in thermal counterflow, above which the 7 
factor suddenly increases, the so-called T1-T2 transition. 
The existence of this transition depends on the details 
of channel cross-section in counterflow (according to the 
data reviewed by Tough two states are observed if the 
aspect ratio is of order one and the smallest dimension 
is less than about 1 mm; a single state is observed if the 
aspect ratio is much larger than 1 or the smallest dimen¬ 
sion is larger than 1 mm). In superflow on the contrary. 
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in channels of different cross-section, only one regime of 
turbulence has been observed. These facts may be re¬ 
lated to the different dynamics of the normal component 
in the two systems. 

Theoretical and computational work has recognized 
that earlier simulations were unrealistic in at least one 
important respect: they related to an unbounded volume 
of helium, and to a situation where the velocity of the 
normal fluid (assumed laminar) and the mean velocity of 
the superfluid are spatially uniform. In practice, thermal 
counterflow takes place almost always in a channel of fi¬ 
nite cross-section, so that, at the very least, the no-slip 
condition for the viscous normal fluid at a solid bound¬ 
ary must lead to spatial non-uniformity in the velocity of 
the normal fluid. Furthermore, the Reynolds number for 
the flowing normal fluid is typically quite large, so that 
it is questionable whether the flow of the normal fluid 
remains laminar. 

The results of more realistic theoretical and computa¬ 
tional work can be summarized as follows. The stability 
of laminar normal-fluid flow has been studied by Melotte 
and Barenghi , although with assumptions that have 
turned out to be unrealistic (the normal-fluid velocity 
profile remains approximately parabolic in the presence 
of mutual friction, and the line density remains unper¬ 
turbed and spatially uniform). The properties of the 
vortex tangle have been studied when the normal fluid 
has a prescribed laminar parabolic profile by Aarts and 
De Waele 1^. Adachi and coauthors a and Baggaley 
and Laurie [^; and when the normal fluid has imposed 
on it a prescribed classical turbulent flow profile by Bag¬ 
galey and Laizet 0 . In the former case the line density 
turns out to be quite inhomogeneous, although the spa¬ 
tially averaged value of 7 is not seriously affected. In the 
latter case the value of 7 is increased, suggesting that the 
transitions at which 7 is observed to increase are associ¬ 
ated with transitions to normal-fluid turbulence. How¬ 
ever, this work is still unrealistic for two reasons: the 
motion in the two fluids, coupled by mutual friction, is 
not treated in a dynamically self-consistent way; and it 
is assumed that there is no pinning of vortex lines at a 
solid boundary. 

We note at this point that our recent experimental 
work on bellows-driven superflow Q showed that the 
steady-state spatially-averaged line density is hardly any 
different from that observed in thermal counterflow at the 
same relative velocity v. This is in spite of the fact that 
the normal fluid velocity, relative to the channel walls, 
is very different in the two cases. This suggests that the 
average line density is insensitive not only to any inhomo¬ 
geneity in the vortex tangle, as suggested by simulations, 
but perhaps also, contrary to the simulations, to any tur¬ 
bulence in the normal fluid, since any turbulence in the 
normal fluid might be expected to be different in the two 
types of flow. We are led to conclude that detailed in¬ 
formation about counterflow turbulence must come from 
experiments other than those that measure average line 
densities in steady-flows: two possible directions are the 


direct visualization of the turbulence and a study of the 
decaying line density after the flow is switched off, which 
we pursue in our work. 

On the visualization side, experiments on ther¬ 
mal counterflow have been reported using as tracers 
both micron-sized hydrogen or deuterium particles and 
metastable He 2 excimer molecules [32[. In the former 
case at relatively low heat currents bimodal distributions 
of vertical velo city have been first measured by Paoletti 
and coworkers [^, indicating that some of the particles 
move in the direction of the normal fluid, while others are 
trapped on vortices and move with the tangle, with a ve¬ 
locity generally different from that of the superfluid. At 
larger heat currents, where effects relevant to our present 
studies might be seen, the particles experience frequent 
trapping and de-trapping: vertical velocity distribution 
changes to a monovalued one [sj and interpretation is in 
general harder. A wealth of statistics of particle velocity 
and acceleration has been produced by the Prague group 
in counterflow in different heat current regimes |^, RO, |35| , 
showing that the character of particle dynamics changes 
from classical to quantum as the length scale investigated 
is reduced from well above to well below the inter-vortex 
separation. Additionally, there has also been some indi¬ 
cation that vortical structures exist on scales larger than 
the intervortex spacing 0- The excimer molecules in¬ 
stead are useful because they track only the normal fluid 
(at temperatures above IK). The use of these molecules 
is still at an early stage of development, but Guo and 
coworkers have already produced evidence that the nor¬ 
mal fluid does become turbulent above a critical veloc¬ 
ity [ 3 ^. Very recently, Guo’s group have been studying 
the time evolution of thin lines of excimer molecules pro¬ 
duced in counterflowing helium, from which they can de¬ 
duce the flow of the normal fluid in greater detail j^ . 
We will further refer to this work in the Discussion ses¬ 
sion. 

We now turn to the main topic of this paper, i.e. the 
decay of vortex line density in superflow and counter- 
flow. The phenomenological Vinen equations describing 
the growth and decay of line density in counterflow turbu¬ 
lence [22[ predict that the free decay should be described 
by the equation 


where X 2 is a dimensionless temperature-de pen dent pa¬ 
rameter proportional to mutual friction [^, , related 

closely to the effective kinematic viscosity v'. It follows 
that 



where Lq is the line density at time t = 0. Recent sim¬ 
ulations of Mineda et al. based on the assumption that 
counterflow turbulence is homogeneous and occurs in an 
unbounded medium with laminar flow of the normal fluid 
have confirmed this prediction [^. However, it is now 
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well-known (2^, liol - E^ , that the experimentally observed 
decay is quite different. Although there is for a very short 
time a rapid decay that might be described by Eq. o, 
the decay then slows down for a time and may actually 
reverse (there is a “bump” in the decay curve), while at 
large times the decay is observed to continue as 
rather than as Note that in the rest of this paper 
we shall use the jargon word “bump” to refer to both the 
case of actual increase in vortex line density during decay 
(i.e. change of sing in dL/dt) and also to the simple slow 
down in decay rate associated to the presence of a point 
of inffexion (i.e. a change of sing in and we will 

be more specific if needed. 

It is now widely accepted that the time-dependence 
as is associated with the decay through a quasi- 

classical Richardson-Kolmogorov cascade of coupled 
(superfiuid-normal fiuid) energy containing eddies, the 
size of which is determined by and limited W the di¬ 
mensions of the containing channel p, [ij, . It has 
been shown that this decay is described in detail by the 
equation 


m = 


D{3Cf/^ 


it - 


( 5 ) 



FIG. 1. (Color online) Flow channel for turbulent superflow 
decay studies, to scale, with dimensions in mm (left). Helium 
flows vertically up, provided by low temperature bellows next 
to the channel. Superleak filters allow flow-through of super- 
fluid component only. Second-sound probes measure the total 
length of quantized vortices per unit volume, Eq. 0. The 
channel can be modified to host thermal counterflow too, as 
illustrated by sketches on the right, where S and N stand for 
superfluid and normal components of He-II. The channel ex¬ 
ists in two variants, with internal square cross-section 7 mm 
and 10 mm in side, coded D7 and DIO in the article. 


where D is taken as the channel width, C ^ 1.5 is the 
classical Kolmogorov-41 constant, and to denotes the vir¬ 
tual time origin. 

Behaviour described by Eq. m was first observe d by 
the Oregon group in the decay of grid turbulence |46l-[4q , 
where the formation of large eddies can be understood in 
classical terms. However, formation of large-scale clas¬ 
sical eddies in the decay of counterfiow turbulence must 
be less straightforward than is the case with grid tur¬ 
bulence. Are they formed out of large scale eddies al¬ 
ready present in the steady state? Or are they gener¬ 
ated from scratch during the early stages of the decay? 
Presumably, the processes involved in this formation are 
reflected in the early stages of the decay of line density, 
but interpretation is hard. There has been some spec¬ 
ulation about these early stages, particularly about the 
origin of the “bump” [^, to which we shall refer later, 
but there is as yet no agreed explanation. Eurther sub¬ 
stantial progress must probably await the results of vi¬ 
sualization experiments, perhaps backed up by more re¬ 
alistic numerical simulations. However, we argue that if 
indeed the large eddies responsible for the decay 

are formed out large-scale eddies present in the steady 
state, then bellows-driven superfiow could exhibit very 
different features in its decay from that seen with ther¬ 
mal counterfiow. We have therefore undertaken a study 
of the decay of such bellow-driven superfiow with a di¬ 
rect comparison of thermal counterfiow studied in the the 
same channel and conditions, and we present the results 
in this paper. As we shall see, the two cases do indeed 
exhibit different forms of decay. 


III. THE EXPERIMENT 
A. Apparatus and method 

A full description of the mechanically driven superfiow 
apparatus and the measurement technique was given in 
Ref. 0 . Here we briefly recall only the essential features 
for convenience, focussing on the aspects more specific 
to the decay studies. A drawing of the flow channel is 
in Eig. [H Two vertical brass flow channels have been 
used in this experiment. The test section is 105 mm in 
length and has an internal square cross-section of side 
7 mm and 10 mm, therefore with a factor 2 change in 
cross-sectional area (coded D7 and DIO in this article). 
The channel ends are plugged by sintered-silver superleak 
disks, each 2 mm thick, 16 mm diameter, with 1/2 filling 
fraction; these discs serve to prevent any net flow of the 
viscous normal component on time scales relevant to the 
experiment. 

The superfiow is driven by a low temperature stainless 
steel bellows immersed in the open cryostat bath and op¬ 
erated through a shaft by a computer-controlled precision 
motor located above the cryostat at room temperature. 
The velocity of the flow in the channel is inferred from a 
measurement of the rate of compression of the calibrated 
bellows volume, assuming that the helium is incompress¬ 
ible. Counterfiow is studied in the traditional way, as 
in previous Prague experiments \^: we used the 
same channel of the superfiow experiment (DIO only), 
installed it horizontally in the cryostat, removed one su- 
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perleak to open the channel to the bath, and fitted the 
other end with a close cap which hosts a 50 heater 
resistor (sketch in Fig. [T]). The dissipated power is mea¬ 
sured continuously, by independent measure of voltage 
and current. 

Turbulence is detected by measuring the extra attenu¬ 
ation of second-sound caused by scattering of normal- 
component thermal excitations by the vortex lines. 
Second-sound is generated and detected by a pair of vi¬ 
brating porous membranes located in the walls of the 
channel at its mid-point (see Fig. [1]); the second-sound 
travels across the channel, which acts as a resonator. The 
time-dependent attenuated amplitude of second-sound at 
resonance a{t) can be related to the instantaneous total 
length of quantized vortex lines per unit volume, L{t) 
through the equation: 


^ GnAfo 
Bn 



(6) 


where ao and A/o are the amplitude and full width at half 
maximum of the second-sound amplitude resonant curve 
for quiescent helium, and B is the mutual friction coeffi¬ 
cient of order unity, tabulated in Ref. (the frequency 
dependence of B can be neglected in this experiment 
since we perform measurements only with a single low 
frequency mode). The attenuation of second-sound mea¬ 
sures the length of vortex line per unit volume weighted 
by a factor sin^ 6>, where 0 is the angle between any el¬ 
ement of vortex line and the direction of propagation of 
the second-sound. To derive Eq. ^ Q, the distribution 
of vortex lines is assumed to be fully random and spa¬ 
tially uniform; the opposite limiting case where the lines 
are instead assumed to be fully polarized, i.e. forming 
vortex rings lying in planes perpendicular to the flow di¬ 
rection, leads to a version of Eq. m a factor 4/3 higher. 
Therefore if the real vortex line distribution is not known, 
the use of Eq. can lead to errors in L{t) of at most 
33%. We will return to how this aspect may affect results 
during discussion. 

The decay process is too fast in time to allow sampling 
of a full resonance curve at any point during the decay. 
However, our studies of the steady-state have shown that 
the second-sound resonant frequency is not significantly 
affected by the extra attenuation, so that it is sufficient to 
set the second-sound frequency on resonance and observe 
this resonant response as a function of time, as shown 
in Eig. [21 This response accurately reflects the chang¬ 
ing attenuation only if the natural response time of the 
second-sound resonator (of order the inverse line-width) 
is sufficiently short; in practice A/“^ is about 10 ms at 
the start of the decay and about 100 ms for quiescent 
helium, values that are small enough for our purposes. 
Amplitude times-series are sampled at 60 Hz, and each 
sample is averaged with a lock-in amplifier with time con¬ 
stant of 10 ms. Temperature control in our cryostat is 
of order 0.1 mK, therefore sufficient to ensure that tem¬ 
perature drift cannot cause significant drift away from 
resonance, as we have experimentally verified. 



FIG. 2. (Color online) Example of a single raw measurement 
of superflow turbulence decay in D7 channel, at T = 1.26 K. 
The second-sound amplitude a{t) (black solid line) is plot¬ 
ted as a function of time. At t = 0 helium is quiescent, at 
t = 30 s a steady-flow of velocity 11.3 cm/s is produced by the 
bellows, and at t = 51.5 s it is suddenly switched off, result¬ 
ing in the following decay. The reference signal ao entering 
the calculation of vortex line density, Eq. is the average 
of the last 20 s of the amplitude signal. The red solid line 
(right axis) shows the bellows displacement. Inset: a zoom 
around the switch-off time, showing the second-sound reacts 
immediately after the halting of the drive. Time resolution of 
our measurements is about 16 ms. 


To create a turbulent steady-state, the bellows are 
compressed at a constant rate for about 20 s (Eig. [2]), 
after which compression is suddenly stopped. The po¬ 
sition of the movable end of the bellows is recorded by 
the motor encoder, which shows that sudden stopping 
is achieved in less than 10 ms (inset). The decay of line 
density is followed for 200 s, after which the amplitude of 
the second-sound signal has reached a statistically steady 
value; a further 30 s is allowed to elapse before a new mea¬ 
surement is made. The second-sound amplitude before 
and long after the steady flow was observed to be gener¬ 
ally slightly different. This effect - related to a varying 
remanent amount of vortex line in the sample - is stud¬ 
ied statistically in Section IHIBI and has some bearing on 
the interpretation of results. It is not clear whether the 
beginning or the end of the amplitude time series ought 
to be used for ao in Eq. we calculate it averaging the 
last 20 s. The measurement protocol for thermal coun¬ 
terflow decays is very similar. 

The parameter space covered by the experiment is 
summarized in Table [H We have performed mechanically- 
driven superflow and thermally driven counterflow decay 
measurements, in the two channels D7 and DIO, in the 
temperature range between 1.25 K to 2.10 K, and for ini¬ 
tial steady-state line densities, Lq, spaced almost exactly 
one decade apart: 10^ cm“^, 10^ cm“^ and 10^ cm“^ (we 
shall refer to the decay curves corresponding to these ini¬ 
tial line densities as 1/6, L5 and 1/4). The table shows 
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Mechanical super flow 

Thermal counterfiow 



D7 

DIO 

DIO 

T[K] 

Lo [cm“^] 

Vs [cm/s] 

Vs [cm/s] 

Vns [cm/s] 


10® 

11.3 

/ 

/ 

1.25 

10® 

3.6 

4.25 

/ 


10^ 

/ 

1.3 

/ 


10® 

7.7 

7.4 

8.45 

1.45 

10® 

2.6 

2.7 

3.0 


10^ 

0.86 

1.0 

1.0 


10® 

5.7 

5.3 

7.1 

1.65 

10® 

1.9 

2.0 

2.5 


10^ 

0.57 

1.0 

0.9 


10® 

4.9 

/ 

/ 

1.75 

10® 

1.7 

/ 

/ 


10^ 

0.6 

/ 

/ 


/ 

/ 

/ 

/ 

2.10 

10® 

0.4 

/ 

1.50 


10^ 

0.15 

/ 

0.51 


TABLE L Overview of the parameter space explored by 
the experiment. Steady-state line density and corresponding 
mean superflow and counterfiow velocity, shown per temper¬ 
ature and channel width (7 mm and 10 mm side of square 
cross-section). 


the corresponding superflow and counterfiow velocity in 
the steady-state, with experimental uncertainty of 3%. 

For every combination of temperature and starting line 
density we have measured typically 150 individual de¬ 
cays, under nominally identical experimental conditions 
and we have ensemble-averaged these samples, by lin¬ 
early interpolating each one onto a 100 Hz time-series 
and averaging point-wise. Decay signals are checked in¬ 
dividually and rare anomalous ones are discarded from 
averaging. The improvement of the averaged signal with 
respect to a single sample is demonstrated in Fig.[3l Av¬ 
eraging over a large ensemble has proven essential for 
our study, allowing us to resolve 6 orders of magnitude 
of decay on L. 


B. The remanent vortex line density 

Our second-sound attenuation technique cannot pro¬ 
vide an absolute measurement of vortex line density, be¬ 
cause the reference second-sound signal ao in Eq. m may 
itself be attenuated by remanent vortex lines persisting in 
quiescent helium after the turbulence has decayed. Vor¬ 
tex lines are indeed expected to survive because they ^n 
effectively to any surface, due to their A-size core [l|. 
The absolute vortex line density in our sample is therefore 
Labs = Lrem+L, whcrc Lrem IS thc nou-measurablc rema¬ 
nent density of vortex lines hidden in ao- However, Lrem 
is expected to be negligibly small compared to L during 
most of the decay. An estimate provided by the work of 



FIG. 3. (Color online). Demonstration of the improvement of 
signal-to-noise ratio by ensemble average of 150 decay sam¬ 
ples. Data relates to pure superflow, D7, T = 1.26 K, L6. 
The noisier signal in solid black line is a single sample ran¬ 
domly chosen from the batch of 150. Time is rescaled here, 
with t = 0 marking the instant when the bellows stops. On 
the y-axis is the vortex line density calculated from the data 
in Fig. [2] via Eq. ©. The average signal only looks noisy at 
late times because it averages to zero and is plotted in loga¬ 
rithmic scale. We can resolve 6 orders of magnitude in decay 
of vortex line density. 



FIG. 4. (Color online) The dependence of the standard de¬ 
viation of the remanent vortex line density Lr (vortex lines 
remaining in the channel after a decay process) with temper¬ 
ature and steady-state line density. Data relates to the D7 
channel, DIO being similar. Lr is calculated with Eq. © 
and its standard deviation is typically from a batch of 150 
decays. The distribution of Lr in a typical batch (here L5 at 
T = 1.65 K) is shows inset, roughly with a Gaussian distri¬ 
bution with zero mean. 










































7 




02468 10 02468 10 02468 10 

t [s] t [s] t [s] 


FIG. 5. (Color online) Decay of vortex line density of pure superflow turbulence. Plots of L (in log scale) versus t for the D7 
channel, at different temperatures. Each panel groups decays from the same steady-state density, at different temperatures. A 
version of this Figure in log-log coordinates is in the Supplemental Material (SM). 


Awschalom and Schwarz [5l| leads to Lrem = 72 cm“^ 
and 36 cni“^ for the D7 and DIO channels, respectively. 
What the experiment can access, however, is the varia¬ 
tion of L in the quiescent state between one decay and 
another, which constitutes a measurement of the rema¬ 
nent vorticity above the unmeasurable I/rem floor. We 
calculate this quantity as 

Lr = -T/gtart -7/end 5 (7) 

where this difference is obtained from the change in aver¬ 
age amplitude of the second-sound signal averaged over a 
period of 20 s in the quiescent states before and after the 
imposition of a flow. We have studied the distribution 
of Lr in the batch of 150 decays across all parameter 
space. An example for D7, T = 1.65 K and L5 is given in 
the inset of Fig. IH showing roughly a Gaussian distribu¬ 
tion centered around zero; this distribution is common to 
all batches. Notice that Lr can be positive or negative, 
and that the extent of its variation cannot be accounted 
for by changes in helium temperature from one decay to 
another (which are too small for this effect); therefore 
we attribute this effect to varying remanent vortex line 
density in the sample. The standard deviation of this 
distribution, SD{Lr), is given in the main plot of Fig. |4] 
as a function of T and Lq for the D7 channel, DIO being 
similar. From this we learn that SD{Lr) is not corre¬ 
lated to I/O and that with the exception of the points 


at T = 2.1 K where measurements are more difficult, 
it varies from about 100 cm“^ to 300 cm“^. We guess 
therefore that an upper limit on the absolute remanent 
vortex line density in the sample is roughly SD (Lr). 

We do not know precisely what effect remanent vor¬ 
tices might have on the rate of decay of a vortex tangle, 
compared to a decay where vortex lines are totally anni¬ 
hilated, a condition which, as we have explained, cannot 
be achieved experimentally due to pinning. However we 
have indication that when the counterflow heat current 
is reduced not to zero, but to a finite subcritical value, 
the decay shape is somewhat affected (Figd) - this sug¬ 
gests that the decay rate may be affected by the amount 
of residual vortices. A detailed study of these effects was 
carried out for thermal counterflow [22| and we plan a 
similar approach for mechanical superflow. 

In addition to these effects a further complication for 
interpretation arises at low densities: for T < 2 K our 
averaged estimated residual density SD{Lr) ^ 250 cm“^ 
corresponds to a ratio of line separation to channel width 
of about 0.1 in D7, at which wall effects might well start 
to be important. For these combined reasons, we think 
it wise to consider our results only up to times for which 
L ^ SD(Li^), despite the fact that our ensemble-average 
curves resolve a longer decay process. For example, for 
the decay in Fig. [3] for which SD{Lr) ^ 100 m“^ we 
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FIG. 6. (Color online) Decay of vortex line density of pure superflow turbulence, as in Fig. [5] but for channel DIO. Features 
are similar to D7, except that L6 curves display a short non-monotonic behavior and the L5 curve at 1.25 K has a change in 
the decay rate with a point of inflexion. This feature, referred to in the text as a “bump”, is an exception in superflow decay 
and the rule in counterflow decay, shown in later figures. 


would limit our consideration up to about t = 90 s. This 
is taken into account in the rest of the article. 


IV. EXPERIMENTAL RESULTS 

An overview of the first 10 s of all our experimentally 
observed superflow decays in the D7 and DIO channels is 
shown in Figs. [5] and [6] respectively. Each panel groups 
decays from the same initial density and different tem¬ 
perature. These plots serve to demonstrate the main 
features of the decays at a glance; details on shorter and 
longer time-ranges will be presented in due course. An 
overview of the full time range in log-log coordinates is 
available for D7 in the Supplemental Material (SM), with 
DIO being similar. Each decay shown here in Eigs.[5]and 
[6] displays an initial fast rate (see also Eig. [H]for details) 
which expires within roughly 0.5 and 2 s depending on 
conditions, followed by a slower process which continues 
for a longer time (see also Eigs.[9]and[TT]for details). Eor 
short time behavior, we notice that qualitatively but sys¬ 
tematically the ratio of Lq to Lx at which the fast process 
changes to a slow process decreases with Lq for fixed T, 
and increases with T for fixed Lq, in both channels. The 
time at which Lx occurs increases systematically with 
decreasing Lq for fixed T and increases with T for fixed 


Lo. 

In D10-L6 only (Eig. [6]- left panel), the fast and slow 
regimes are joined by an intermediate one, with an in¬ 
version of the decay rate, i.e. a change of sing of dLjdt. 
In D10-L5 at T = 1.25 K, instead, we observe only a 
slow down of the decay rate without an increase in L, 
but with the presence of a point of inflexion across which 
changes sign. We shall refer henceforth to these 
features as a “bump” in the decay curve. Eor simplic¬ 
ity we shall call “bump” both the increase in L during 
the decay and the presence of a point of inflexion, and 
we shall be more specific when needed. The bump, ob¬ 
served only in these few circumstances in superflow, is 
instead seen always in counterflow in the stricter sense 
of non-monotonic behaviour, as demonstrated in Eig. [7| 
and related Supplemental Material (SM). The existence 
of a non-monotonic bump in counterfiow was already de¬ 
tected in earlier experiments j42l-[4^ and here, within the 
range of investigated parameters, we confirm it. There 
were however instances in past experiments, including 
the once just cited, when a point of inflexion was ob¬ 
served instead of non-monotonic behaviour, as in the first 
experiments on the decaying thermal counterfiow by Vi- 
nen or more recently in the Prague decays [U , 
for decays originating from steady-state counterfiow gen¬ 
erated by larger heat fluxes than reached here. What 
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Fig. [ 7 ] adds to previous studies is high statistics for coun- 
terfiow (ensemble average of 150 decays instead of single 
curve, i.e. the same standards as for our superflow) and, 
especially the fact that we can compare superflow and 
counterflow strictly under the same conditions: channel, 
temperature, and initial vortex line density. 

Fig. [ 7 | shows that the bump in counterflow is broad¬ 
ened with decreasing Lq, and that in superflow it occurs 
only for L6, although the details of the decay remain dif¬ 
ferent in the two cases. The lighter curves in the left 
panel are decays when the heat current is reduced, not 
to zero, but to a small value (10 mW) below the criti¬ 
cal value for transition to quantum turbulence. The idea 
here was to observe if leaving a subcritical heat current 
in the channel may help in reducing the level of residual 
vortex lines, “washing them away”, in the spirit of the 
systematic study in Ref. [22|. Whilst we did not observe 
a change in the residual second-sound attenuation at the 
end of the decay process, we did measure differences dur¬ 
ing the rest of the decay, as shown. At late times the line 
density is lower in case of non-zero heat current, and the 
“bump” is somewhat reduced. 

Let us now focus on the fast initial decay regime ob¬ 
served at short times, emphasized in Fig. [8] for counter- 
flow and superflow in DIO at T = 1.65 K. In the fig¬ 
ure we compare this decay regime with the prediction 
of the Vinen model in Eq. (|1]) relating to the decay of 
a fully random tangle unbounded by walls, by recasting 
the equation as follows: 



where the quantity t<iX 2 /27r is given in Table [TTl We no¬ 
tice that the Vinen decay is not generally observed, and 
when it exists, is followed for only a very short time, at 
most 1 s. The pure superflow is more Vinen-like than the 
counterflow, and departures increase with increasing Lq 
in both cases. The situation is similar at other temper¬ 
atures, with departure from Vinen behavior increasing 
with decreasing temperature and increasing initial den¬ 
sity (data at T = 1.45 K are provided as SM). In Fig. [8] 
we notice also that the decay does not start abruptly as 
predicted by Eq. but there is some rounding imme¬ 
diately after t = 0. This rounding is more pronounced in 
superflow than in counterflow, lasting at most some 200 
ms, with the tendency to increase with initial line den¬ 
sity. This latter fact leads us to favour the explanation 
that the rounding may be an instrumental effect. The 
switching off of counterflow is controlled fully electron¬ 
ically, whilst the bellows is a mechanical system which 
may introduce secondary lagging effects. Although the 
bellows motor encoder at room temperature does show 
that the system comes to rest to within 10 ms (Eig. [2]), 
the actual flow may not stop abruptly, for reasons such 
as finite compressibility, the finite time for expiration of 
pressure gradients (which would increase with bellows ve¬ 
locity as observed), the expiration of thermal gradients 
occurring because of the presence of superleaks causing 


T[K] 

[cm^s 

1.25 

2.2 X 10“^ 

1.45 

3.1 X 10“^ 

1.65 

4.5 X 10“^ 

1.75 

5.3 X 10“^ 

2.10 

9.3 X 10“^ 


TABLE 11. Values of obtained averaging the values 

computed by Schwarz and those given by the theory of 
Vinen and Niemela Q]- This quantity occurs in Eq. which 
is plotted in Fig. [8l 

an increase of temperature in the bellows and a decrease 
in the channel due to change in superfluid/normal density 
ratio (see Ref. Q for more detailed discussion), and so on. 
In principle however one should not exclude the possibil¬ 
ity that the flow actually stops at t = 0 and therefore the 
rounding would be explained by some incompleteness of 
the model in Eq. (|1]). At any rate, our time resolution of 
^ 16 ms, limited by the intrinsic physics of second-sound 
resonance, does not allow us to study this physical pro¬ 
cess in greater depth. We therefore concentrate on the 
decay process after the first, say, 100 ms have elapsed. 

The situation at late times is summarized in Eig. O 
for the case of superflow in D7, with starting density L5 
at different temperatures (DIO version in SM). Although 
we know from section ME that we should handle with 
care the decay process when L becomes comparable with 
or smaller than its mean remanent value, in this figure 
we demonstrate that a single power law of the form L oc 

— represents the data from the first few seconds 

to the end of the range. This is the behavior predicted 
by the quasi-classical model in Eq. 

Since Eq. (|5]) has a virtual time-origin to, this must be 
subtracted from the time axis in a logarithmic plot. To 
estimate it we recast Eq. as follows 

and obtain to as the intercept of the linear part of the 
plot with the time axis, as demonstrated in the inset 
of Eig. [9] by dashed lines. The time at which the t“^/^ 
behavior onsets, i.e. the “saturation time” tg (“satura¬ 
tion” refers to the attainment of the condition of large 
eddies reaching their maximum size limited by the chan¬ 
nel width D, required in deriving Eq. ®), is plotted for 
different flows in Eig. [10] against the Reynolds number, 
defined with the mean velocity u, the channel width D, 
and the kinematic viscosity of the normal component 
following Ref. [i^- Allowing for some scatter, all exper¬ 
iments produce roughly a scaling tg oc Re~^. 

We continue the analysis of the late time behavior by 
contrasting superflow and counterflow in Eig. [TTl foi* fhe 
same temperature, T = 1.65 K and channel, DIO, observ¬ 
ing that in both cases there is linear range as predicted by 
Eq. ([9]) extending as far as 80 s- the similar T = 1.45 K 
case is in SM. This equation also predicts that (i) curves 








10 




0 5 10 15 20 25 


t[s] 


FIG. 7. (Color online) Comparison of decays observed for thermal counterflow with those observed for bellows-driven superflow 
for different initial line densities, at the same temperature, T = 1.65 K, and in the same channel, DIO. Line densities are 
normalized to the initial line densities. For counterflow, lighter lines correspond to decays when the heat current was reduced 
to a sub-critical value of 10 mW. The temporary increase in vortex line density during decay (the “bump”) is observed 
systematically in thermal counterflow, but only in the DIO channel and at high density in mechanical superflow. Similar plots 
at T = 1.45 K are in the Supplemental Material (SM), including also the thermal counterflow decay at T = 2.10 K which does 
not have a counterpart in superflow for the same channel size. 
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FIG. 8. (Color online) Comparison of decays observed for thermal counterflow with those observed for bellows-driven superflow 
for different initial line densities at very small times for T = 1.65 K. Plots of 1/L — I/Lq against t according to Eq. (jS]). The 
Vinen decay of a random tangle is observed for a very short transient and is followed more closely by mechanical superflow 
than thermal counterflow, and, in both cases, by decays originating from lower densities. (The T = 1.45 K data are in SM) 
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FIG. 9. (Color online) Superflow decay in the D7 channel 
from initial line density Lo = 10^ cm“^ at different temper¬ 
atures. In the main plot L and t are plotted in logarithmic 
axis, demonstrating that superflow decays with the power law 
L (X (t — (dashed line) for a very large fraction of 

time, except for the first few seconds, according to the quasi- 
classical single fluid decay model in Eq. (|5)). The x-axis re¬ 
quires the subtraction of the virtual time origin to obtained 
from the inset, as described in the text. Similar plots from 
different initial density and for the DIO channel are in SM 
and confirm the same conclusions. 



Re = vD/v 


FIG. 10. (Color online) The time for the onset of L oc 
behavior (saturation time, see text) as a function of the flow’s 
Reynolds number. The present experiments are compared 
with an older Prague thermal counterflow experiment in a 
round channel in and with the Oregon towed grid exper¬ 
iments (the dashed-line is an interpolation of their experi¬ 
mental points) which measures the turbulence decay behind 
a moving grid [4^ . 


from different initial densities should collapse: our data 
confirms that decays from different Lq have very similar 
slope in the regime, but the collapse is observed 

only in counterflow and not in superflow, (ii) The slope 
should scale as despite the fact that we have only 

two channel sizes, we checked this prediction and found 
that experimental slopes scale by a substantially larger 
extent than expected. In these two respects therefore our 
results differ from the prediction. 

Fitting Eq. (|9|) to the linear part of the curves in Fig.fTT] 
allows us to extract the effective viscosity as is custom¬ 
arily done [illSil. In Fig. [121 together with the data 
from the decay of turbulence past the towed grid [52| , we 
thus plot iy'{T) for all our flows. The effective viscosity 
measured from decays starting from different vortex line 
density Lq is found to slightly vary, but in a manner un¬ 
correlated with Lq. In any case, we do not expect u' to be 
dependent on Lq because we understand to be a robust 
property of the flow independent of flow details, as it is 
shown by consistent values coming from different decay 
experiments [2^ and also from steady-flows, hence from 
an entirely independent approach [53|. For these reasons, 
we think it justified to average z/' originating from differ¬ 
ent I/O - which is reflected in the error bars in Fig.[T2j We 
notice that, despite the data from the three experiments 
at the same temperature being roughly consistent within 
error bars, the superflow data for DIO lie systematically 
lower than those for D7, which puts pressure on the un¬ 
derstanding that ly' should not depend on flow details. 
At any rate, as we have discussed in our works which 
include determination of effective viscosity [bsl, [bj , the 
absolute value of u' obtained from decay measurements 
can have uncertainties up to a factor of 4, depending on 
how strictly the assumption on saturation of the large 
eddies size attains. Accurate experimental determina¬ 
tion of z/' to better than a factor of 4 remains therefore 
still a challenge. 


V. DISCUSSION 

All our experimentally observed decays, both in su¬ 
perflow and counterflow, have an initial regime when the 
line density drops rapidly (for sufficiently low density it 
approaches for a short time the Vinen law L (x es¬ 
pecially in superflow), and a final slower regime of the 
form L oc The first, fast regime, lasts at most 1 

s (Fig. [8]) and is responsible for the loss of a large frac¬ 
tion of vortex lines, even up to about 99% (see, e.g., the 
T = 1.75 K curve in Fig. |5l left panel). The fact that 
the decay rate increases with increasing temperature for 
fixed initial density may be due to rising value of the 
prefactor /^X2/27r in Eq. (|3j) with temperature, as shown 
in Table ini which in turn must be related to the increase 
in mutual friction. The second, slow regime, lasts for the 
majority of the decay process and causes the loss of the 
remaining vortex lines (except the remanent ones). 

We have compared our observations with two available 
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FIG. 11. (Color online). Comparison of decays observed for thermal counterflow with those observed for bellows-driven 
superflow for different initial line densities at T = 1.65 K in DIO. Plots of against t as according to Eq.(|9}. Linear 

behavior corresponds to the quasi-classical L oc observed in both flows. The T = 1.45 K version is in SM. 
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FIG. 12. (Color online) The effective viscosity of the “quasi- 
classical” model (Eq. (O) as a function of temperature. Eor 
our three experiments, at each given temperature, u' has been 
averaged from decays of different initial density (L4, L5 and 
L6) , resulting in the shown error bars. Eor comparison we 
add also data from the Oregon towed grid experiment [H^] • 


analytic models: (i) the Vinen equation for a random tan¬ 
gle - and found it to apply in a limited regime and (ii) the 
quasi-classical model for the decay of large eddies - and 
found it to apply more generally although not fulfilling 
all predictions. 

Our results generally confirm the understanding 0 
Issl 0 that a general quantum turbulence tangle con¬ 
sists of an random system of vortex lines which decays 


fast by mutual annihilation of lines, a fraction of which 
is organized in bundles giving rise rise to eddies of var¬ 
ious sizes, up to sizes comparable to the channel width, 
which decay more slowly. Alternatively, we can imagine 
there being a more spatially homogeneous regime with 
an energy spectrum that is Kolmogorov in form except 
for an initially enhanced energy at wave numbers close 
to the inverse vortex line spacing; i.e a situation in which 
the density of vortex lines is initially larger than is nec¬ 
essary for the dissipation of energy, given by Eq. o , at 
the rate required to match the flow of energy down the 
Kolmogorov cascade. Our experiment cannot establish 
unequivocally whether these eddies exist already in the 
steady-state, although evidence that they do exist comes 
from visualization of counterfiow turbulence both from 
tracing solid particles [T^ and imaging the normal com¬ 
ponent (si, |33- 

What is not predicted by the existing analytic models 
is the occurrence of the “bump”, i.e., the change of sign 
in dLjdt or ^Ljdt^ between the fast and slow regimes, 
which we have observed always in thermal counterflow 
and as an exception in superflow. The second-sound at¬ 
tenuation can become increased during decay for one, or 
for a combination of two reasons: (i) the line length stays 
approximately constant but it is spatially rearranged so 
that the second-sound “sees” a greater fraction (recall 
that second-sound is attenuated only by the projection 
of vortex lines onto a plane perpendicular to second- 
sound’s propagation direction), or (ii) the spatial orien¬ 
tation stays fixed but the line length increases. 

The possibility of(i) has been confirmed numerically 
by Barenghi et al EsI for spatially rearranged ran- 
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dom tangles, but although the simulation gave a qual¬ 
itative result in agreement with observations, the vor¬ 
tex line density in the simulation was about an order of 
magnitude below the lowest available experimentally. We 
note nevertheless that the observed height of the bump 
is never greater than can be accounted for by this mech¬ 
anism (33%). On the other hand, this effect was not ob¬ 
served in the simulations of Mineda and coworkers [39| , 
based on more realistic line densities. However, neither of 
these simulations takes account of the possibility that the 
vortex tangles with which we are dealing are polarised in 
such a way that large scale eddy motion is superimposed 
on the random tangle. 

Option (ii) is also physically possible, and to build 
an argument for it we note that in our measurements 
the bump occurs always in thermal counterflow and only 
at high steady-state velocities in superflow. This sug¬ 
gests that the dynamical state of the normal component 
(laminar, unstable or turbulent) might be relevant to the 
existence or not of the bump, since the average veloc¬ 
ity of the normal component relative to the walls varies 
significantly in the two systems: it can be rather large 
in the case of thermal counterflow, but it is nominally 
zero in the case of superflow. We note however that, 
as suggested in Ref. Q, from considerations on the scal¬ 
ing of critical velocity with channel width from different 
experiments, the normal component is probably not at 
rest in superflow in a large channel, but is set in mo¬ 
tion by mutual friction, the spatially-averaged velocity 
remaining zero. Nevertheless this motion may be rela¬ 
tively slow in comparison with that in thermal counter- 
flow, except at large superflow velocities. But just how 
would the motion of the normal component cause the 
decay inversion? The recent results from the Tallahassee 
group [33 where the normal component is tracked by the 
excimer molecules, give evidence that the normal compo¬ 
nent at a sufficiently high heat current becomes turbu¬ 
lent (specifically at T = 1.83 K for a heat current above 
qc = 80 mW/cm^ at which Lq ^ 10^ cm“^). The sec¬ 
ond order structure functions calculated from the normal 
fluid turbulent velocity fluctuations extend over a fairly 
large range of lengths scales, up to a sizeable fraction of 
the channel size. The structure function can be related 
to a turbulent energy spectrum of the form E{k) oc k~‘^ 
in the steady-state (where k is the wavenumber). During 
the decay the spectrum changes gradually, within about 
3 seconds, into the Kolmogorov form E{k) oc It 

seems possible that the bump is associated with this 
evolution. We are currently exploring this idea with a 
model for the temporal evolution of the energy spectrum, 
and this work will be reported in a future publication. 
Further experimental evidence that the bump is associ¬ 
ated to transition to turbulence in the normal compo¬ 
nent would come by testing thermal counterflow at heat 
currents below the transition. In our system this would 
mean to study decays from initial density sufficiently be¬ 
low I/O = 10^ cm“^, which constitutes a direction for 
future work. 


VI. CONCLUSIONS 

We have presented a comprehensive picture of the tem¬ 
poral decay of vortex line density in quantum turbulence 
produced by mechanically driven superflow through two 
square channels of 7 and 10 mm side. We have cov¬ 
ered a broad parameter space in temperature (1.25 < 
T < 2.10 K) and in steady-state vortex line density 
(10"^ < L < 10^ cm“^). Additionally, at T = 1.45 K 
and T = 1.65 K and for all the same initial densities, we 
have provided, for the first time, a direct comparison of 
mechanical superflow and thermal counterflow, the lat¬ 
ter performed in the same 10 mm wide channel used for 
superflow and under exactly the same experimental con¬ 
ditions. This, together with enhanced accuracy achieved 
thanks to ensemble averaging of up to 150 individual de¬ 
cays placed us in a strong position to compare these flows. 

In an unbound system superflow and counterflow ought 
to display identical physics since they are related by 
Galilean invariance. In practice, the presence of chan¬ 
nel walls will change the physics, at least for the normal 
component which must acquire a profile due to viscous 
drag with the walls. This has been indeed confirmed ex¬ 
perimentally by visualizing the normal fluid flow profile 
using helium excimer molecules, observed to turn from 
laminar to turbulent as heat current increases 371. Addi¬ 
tionally, we know from numerical simulations 3lf, [sl, 
that the normal fluid prohle induces inhomogeneity in 
the distribution of vortex lines across the channel width, 
with the density being enhanced in case of turbulent 
normal fluid prohle. Our experiments with steady-state 
superhow have shown that when comparing with other 
counterhow experiments with channels of different size Q 
and also when measuring in our own superhow channel 
adapted for counterhow [13 , the line density for a given 
relative velocity is essentially insensitive to the normal 
component net how. This is not so for the decay. 

The decay of these two hows is similar in that both 
display an initial fast process where most of the line den¬ 
sity is lost which for sufficiently low density has the form 
L oc and a subsequent slow L oc process where 
the rest of the tangle decays. These two processes have 
been associated respectively with the decay of the ran¬ 
domized and polarized components of the tangle. The 
key diherence however is that in counterhow we invari¬ 
ably observed an inversion of the decay rate between the 
two regimes, which we observed only at high steady-state 
velocities in the wider channel in superhow. This fact, 
hrmly established by high-statistics measurements indi¬ 
cates that the dynamical state of the normal component 
in the steady-state (no net how through the channel in 
superhow and turbulent pipe how in counterhow) has 
consequences for the decay of the tangle. We have spec¬ 
ulated about the reasons, but further work is required if 
these reasons are to command conhdence. We hope that 
this comprehensive set of experimental data describing 
the vortex tangle decay in superhow and counterhow will 
stimulate the development of a still missing full theory of 
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counterflow turbulence. 
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